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FUNDAMENTALS  OF  ION  MOTION  IN 


ELECTRIC  RADIOFREQUENCY  MULTIPOLE  FIELDS 

I.  INTRODUCTION 

The  strong-focusing  effect,  which  utilizes  alternating  electromagnetic  fields,  is  the  basis 
for  high-energy  accelerator  physics1  and  quadrupole  mass  spectrometry.2  Ion  motion 
through  such  electromagnetic  fields,  induced  by  arbitrary-shaped  pole  pieces,3  cannot  be 
described  in  closed  mathematical  form.  However,  the  differential  equations  describing 
ion  motion  through  the  electromagnetic  field  induced  by  voltages  on  an  electrode  struc¬ 
ture  employing  hyperbolic-shaped  electrodes  can  be  written  in  an  analytic  form  which  is 
known  as  the  Mathieu  equation,  the  theory  and  properties  of  which  are  well  known.4'7 
Therefore,  it  is  not  surprising  that  many  mass  filters  are  based  on  electrodes  with  a  hyper¬ 
bolic  geometry. 

The  methods  for  obtaining  the  potential  functions  and  the  differential  equations  of  ion 
motion  in  time-varying  electromagnetic  fields  are  probably  familiar  to  accelerator 
physicists  and  electrical  engineers.  Surprisingly,  a  search  of  the  literature  has  failed  to 
reveal  a  fundamental  and  practical  discussion  of  this  subject  suitable  to  those 
researchers,  such  as  mass  spectrometrists  and  ion  physicists,  who  might  benefit  most  by 
such  a  discussion. 

Research  on  mass  filters  has  been  confined  primarily  to  the  quadrupole.  the  monopole, 
the  three-dimensional  quadrupole,  and  approximations  to  these  three.7'9  It  may  be  possi¬ 
ble  that  electric  radiofrequency  mass  filters  with  desirable  properties  could  be  im¬ 
plemented  if  electrode  geometries  other  than  the  above  were  considered. 


*  E.  D.  Courant,  M.  S.  Livingston,  and  H.  S.  Snyder,  Phys.  Rev.,  88,  1190  (1952);  J.  D.  Blewett,  Phys.  Rev.,  88, 
1197(1952). 

^  W.  Paul  and  H.  Steinwedel,  Z.  Naturforsch.  Teil  A.,  8,  448  (1953);  W.  Paul,  H.  P.  Reinhard,  and  U.  von  Zahn, 
Z.  Phys.  152,  143  (1958). 

3 

The  terms  conductor,  electrodes,  and  pole  pieces  will  be  used  interchangeably. 

*  F.  M.  Arscott,  Periodic  Differential  Equations,  The  Macmillan  Co.,  New  York  (1964). 

®  N.  W.  McLachlan,  Theory  and  Application  of  Mathieu  Functions,  Oxford  University  Press,  New  York  (1947). 

"  Computation  Laboratory,  Tables  Relating  to  Mathieu  Functions.  U.S.  Bureau  of  Standards,  Columbia  University 
Press,  New  York  (1951). 

^  P.  H.  Dawson  (Editor),  Quadrupole  Mass  Spectrometry  and  Its  Applications,  Elsevier,  Amsterdam  (1976). 

®  J.  F.  J.  Todd  and  G.  Lawson,  MTP  International  Review  of  Science.  Physical  Chemistry,  Mass  Spectrometry. 
Series  Two,  Vol.  5,  edited  by  A.  Maccoll,  Butterworths,  London  (1975). 

Q 

J.  E.  Campana,  Int.  J.  Mass  Spectrom.  Ion  Phys.,  S3,  101  (1980)  and  references  cited  within. 
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Previous  work  with  the  quadrupole  mass  filter  electrode  geometry  has  been  directed 
towards  making  circular  cross  section  electrode  approximations  to  “ideal”  hyperbolic 
cross  section  electrodes.  Denison10  has  pointed  out  a  propagation  of  errors  in  the  literature 
that  deal  with  the  correct  positioning  of  cylindrical  electrodes  to  approximate  a  hyper¬ 
bolic  electrostatic  field.  Denison  presented  the  equations  of  ion  motion  in  such  fields  and 
he  reported  variations  in  the  stability  diagrams  resulting  from  the  relative  positioning  of 
the  cylindrical  electrodes  approximating  ideal  hyperbolic  cross  section  electrodes. 
Denison  discusses  the  power  series  expansion  expression  for  the  quadrupole  potential 
distribution  which  is  obtained  by  considering  a  point  charge  located  in  a  quadrupole 
field  formed  by  electrodes  of  line  charge.11  Successive  terms  in  the  series  expression  repre¬ 
sent  quadrupole,  dodecapole,  icosapole  geometries,  etc.  A  commercial  instrument,12 
employing  correcting  electrodes  to  account  for  the  dodecapole  field  distortions  was 
available  for  a  number  of  years  but  ultimately  was  replaced  in  favor  of  electrodes  with 
hyperbolic  cross  sections. 

Although  the  quadrupole  geometry  has  been  used  for  the  mass  filter  in  practice,  there  is 
the  possibility  that  other  geometries  might  give  equivalent  or  superior  resolution  for  a 
given  transmission  and  might  provide  other  advantages  as  well.  New  geometries  could  be 
investigated  empirically,  but  for  practical  reasons  this  task  requires  a  means  of 
evaluating  electrode  geometry  designs  theoretically.  In  this  report,  several  alternatives  to 
the  quadrupole  geometry  are  considered  and  equations  of  ion  motion  in  closed  analytical 
form  are  derived  for  each  of  them.  By  detailed  consideration  of  these  particular 
geometries  it  may  be  possible  that  insight  into  the  characteristics  of  a  useful  mass  filter 
geometry  will  be  gained. 

The  first  objective  of  this  report  is  to  present  the  fundamental  laws  of  physics  and  the 
approximations  which  govern  ion  motion  through  electromagnetic  fields  in  the  electric 
RF  mass  filter.  A  second  objective  is  to  present  a  method  for  finding  the  potential  func¬ 
tion  of  an  arbitrary  “ideal  multipole  field”  and,  subsequently,  the  differential  equations 
of  motion  for  an  ion  traversing  such  a  field.  A  summary  and  discussion  of  forces  and  ef¬ 
fects  which  are  typically  neglected  in  the  physics  of  ion  motion,  but  which  may  be  impor¬ 
tant  under  certain  circumstances,  are  included  for  completeness.  Finally,  the  potential 
functions  and  differential  equations  which  describe  ion  motion  in  the  electric  RF  dipole, 
quadrupole.  hcxapolc.  octapole.  and  decapole  fields  are  derived  using  mathematical 
methods.  This  is  a  first  step  toward  investigating  the  properties  of  alternative  multipole 
geometries  for  mass  filters. 


D.  R.  Denison,/  Vac.  Sci.  Technol.,  8,  266  (1971). 

11  See  H.  Matsuda  and  T.  Matsuo,  Int.  J.  Mass  Speclrom.  Ion  Physics,  24,  107  (1977),  for  a  partial  derivation 
of  the  power  series  expansion. 

*  ^  Hewlett-Packard  Corporation,  Palo  Alto.CA. 


2 


II.  FUNDAMENTAL  PHYSICS  OF  IONS  IN  ELECTROMAGNETIC  FIELDS 

The  electric  RF  strong-focusing  effect  has  been  employed  in  quadrupole  mass  spec¬ 
trometry  for  over  twenty  years.2  9  The  theory  of  operation  for  this  versatile  device  has  been 
discussed  by  many  authors7*9  and  in  all  discussions  it  appears  that  the  fundamental 
physics  underlying  ion  motion  in  electromagnetic  fields  has  been  assumed. 

In  this  section,  a  method  for  finding  the  differential  equation  governing  the  motion  of 
an  ion  through  electromagnetic  fields  with  any  number  of  arbitrary-shaped  pole  pieces  is 
derived.  This  derivation,  while  general,  is  directed  specifically  to  ion  motion  through  an 
electric  RF  field  which  finds  popular  use  as  a  mass  filtering  device. 

The  differential  equation  which  describes  the  motion  of  an  ion  through  a  field  is  ob¬ 
tained  from  Newton’s  law  of  motion 


F  = 


(1) 


where  F  is  the  force  on  the  ion,  tn  is  the  mass  of  the  ion  and  r  is  the  radius  vector  from  the 
origin.  The  force  on  an  ion  with  charge  e  moving  through  an  electric  field  E  and  a 
magnetic  field  B  with  a  velocity  v  is  given  by  the  Lorentz  force  law 

F  =  e(E  +  vxB).  (2) 

It  is  necessary  to  compute  the  E  and  B  fields  from  the  time-varying  potential  imposed  on 
the  pole  pieces  to  find  the  motion  of  the  ions  from  Equations  1  and  2.  These  are  com¬ 
puted  from  the  general  Maxwell  equations: 


V-D=p  VxE  =  -■££. 

dt 

V*B=0  VxH  =  J  +  i£  , 

3t 


(MI) 


where  p  is  the  real  charge  density  in  the  medium  and  J  is  the  current  per  unit  area  in  the 
medium. 

7  W.  Paul  and  H.  Steinwedel,  Z.  Saturforsch.  Teil  A.,  8,  448  (1953);  W.  Paul,  H.  P.  Reinhard,  and  U.  von  Zahn, 
Z  Phys  152.  143(1958). 

7  P.  H.  Dawson  (Editor),  Quadrupole  Mass  Spectrometry  and  Its  Applications,  Elsevier,  Amsterdam  (1976). 

fi 

)  I .  J.  Todd  and  (i.  Lawson,  MTP  International  Review  of  Science.  Physical  Chemistry.  Mass  Spectrometry, 
Series  Two,  Vol.  5,  edited  by  A.  Maccoll,  Butterworths,  London  (1975). 

a 

J.  K.  Cam  pan  a.  Int.  J  Mass  Spectrom.  Ion  Phys.,  33.  101  (1980)  and  references  cited  within. 
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The  del  operator  V  which  appears  in  Equations  Ml  can  be  defined  by 

V<0*d r  =  d<j),  (3) 

where  <0  is  a  scalar  function  of  position,  d0  is  the  total  differential  of  0,  and  dr  is  the 
total  differential  of  the  radius  vector  r.  This  definition  can  be  used  to  find  an  explicit 
representation  of  the  operator  V  in  rectangular  coordinates.  From  elementary  calculus: 


d <t>  =  dx  +  dy  +  ■!!£.  dz. 

dx  dy  9z 


In  rectangular  coordinates  the  displacement  dr  is 

dr  =  idx  +  jdy  +  fidz. 

Since  dx,  dy,  and  dz  are  independent,  in  rectangular  coordinates  V<f>  is  identified  as 


V0 

dx  dy  dz 


and  the  operator  V  is  identified  as 

v=  'LL  , 

9x  9y  3z 


(4) 


where  i,  j.  and  k  are  unit  vectors  along  the  x,  y,  and  z  axis.  Expressions  for  the  Voperator 
in  co-ordinate  systems  other  than  the  rectangular  one  are  given  elsewhere.”  14 

Equations  1  and  2  and  Maxwell’s  equations  Ml  are  fundamental  equations  which 
describe  the  motion  of  an  ion  through  any  electromagnetic  field.  What  do  these  equations 
mean?  Suppose  that  the  charge  density  p(x,y,z,t)  and  the  current  density  J(x.y,z.t)  are 
known  throughout  space  and  for  all  time.  With  V defined  in  Equation  4,  the  divergence, 
curl,  and  the  time  derivative  of  all  the  quantities  in  Equations  Ml  can  be  readily  com¬ 
puted  providing  E(x,y,z,t),  D(x,y,z,t),  B(x,y,z,t),  and  H(x,y,z,t)  are  specified.  Maxwell’s 
equations  assert  that  the  E,  D,  B.and  H  fields  cannot  be  chosen  arbitrarily  and  must  be 
chosen  in  such  a  way  that  they  satisfy  Equations  Ml.  Furthermore,  the  E  and  D  fields  are 
related  and  the  B  and  H  fields  are  related  by  the  properties  of  matter. 

*■*  E.  M.  Pugh  and  E.  W.  Pugh,  Principles  of  Electricity  and  Magnetism,  Addison-Wesley  Publishing  Co.,  Reading, 
MA  <  1960). 

14 

S.  Ramo,  J.  R.  Whinnery,  and  T.  Van  Duzer,  Fields  and  Waves  in  Communication  Electronics,  John  Wiley  and 
Sons,  In-  ;w  York  (1965). 
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The  vector  D  (called  the  electric  displacement  vector)  is  related  to  the  vector  E  by13  14 

D  =  e  E  +  P,  (5) 

where  €o  is  a  constant  called  the  permittivity  of  free  space  and  P  is  a  polarization  vector 
characteristic  of  a  particular  material.  In  linear,  homogeneous,  isotropic  materials,13  14 
P  =  eoxpE(Xe  is  a  dimensionless  quantity  called  the  polarizability  or  electric  susceptibili¬ 
ty)  so  that  Equation  5  becomes 

D  =  (f  >  +  e((X(.)E  =  eE.  <6> 

The  constant  e  in  the  last  equation,  the  permittivity,  depends  only  on  the  medium  in 
which  the  field  exists.  Similarly,  H  (called  the  magnetic  field  intensity  or  magnetic  field 
strength  vector)  is  related  to  B  by13 

B  =  M,H+m,  (7) 

where  fi  is  a  constant  (called  the  permeability  of  free  space)  and  HJ  is  called  a  magnetic 
polarization  vector.15  In  linear,  homogeneous,  isotropic  materials  m  =  fye-H  (the 
dimensionless  quantity  Xm  is  called  the  magnetic  susceptibility  of  the  material  under  con¬ 
sideration)  so  that  Equation  7  becomes 

B  =  /iH.  (8) 

The  constant  n  in  Equation  8,  called  the  permeability,  depends  only  on  the  medium  in 
which  the  field  exists.  Thus  for  linear,  homogeneous,  isotropic  media.  Maxwell's  Equa¬ 
tions  Ml  must  be  satisfied  subject  to  the  constraints  of  Equations  6  and  8. 

With  these  assumptions  (Equations  6  and  8)  about  the  properties  of  the  material  from 
which  pole  pieces  are  constructed  and  for  the  medium  in  which  the  ion  moves.  Maxwell's 
equations  Ml  become:16 

'  ^  E.  M.  Pugh  and  E.  W.  thigh,  Principles  of  Electricity  and  Magnetism,  Addison-Weslcy  Publishing  Co.,  Reading. 
MA  (1960). 

S.  Ramo,  J.  R.  Whinnery,  and  T.  Van  Duzer,  Fields  and  Waves  in  Communication  Electronics,  John  Wiley  and 
Sons,  Inc.,  New  York  (1965). 

Some  authors  write  Equation  7  in  the  form  B  =  y(H+M)  and  call  M  the  magnetization.  This  allows  the  units 
of  H  and  M  to  be  the  same.  The  form  of  Equation  7  is  used  here  to  emphasize  the  similarity  with  Equation  5. 

Here  Maxwell's  equations  are  expressed  in  terms  of  E  and  B  rather  than  D  and  H  since  E  and  B  are  the  vectors 
in  Equation  2. 
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V  •  (eE)  =  p  VxE  =  - 

3t 


V-B  =  0 


J  +  3(eE)_ 
3t 


(M2) 


Since  the  ion  of  interest  is  moving  in  free  space  (vacuum),  e  and  fx  have  the  constant 
values  eo  and  po  in  Equations  6  and  8.  If  the  charges  of  ions  and  the  associated  current 
densities  can  be  neglected  for  the  purpose  of  computing  the  E  and  B  fields  acting  on 
them.  Maxwell’s  Equations  M2  become: 


V'E=0 
V-B=0 

Now  suppose  that  the  E  and  B  fields  are  not  changing  with  time.  For  this  situation.  Max¬ 
well’s  Equations  M3  become: 

Vx  F  =>  0 

V-B  =  0  Vx  B  =  0.  (M4) 


VxE  =  - 


3B 

3t 


VxJL  =  ^e°E* 

u  3t 


(M3) 


The  validity  of  the  different  forms  of  Maxwell’s  Equations  (Ml,  M2,  M3.  M4)  merits 
discussion.  Maxwell’s  Equations  in  the  form  of  Equations  Ml  are  universally  valid  in 
any  medium.  Maxwell’s  Equations  M2  are  valid  without  approximation  for  elec¬ 
tromagnetic  fields  within  a  vacuum  where  the  E  and  D  fields  are  parallel  and  the  B  and 
H  fields  are  parallel.  This  condition  is  fulfilled  when  the  fields  exist  in  linear,  . 
homogeneous,  isotropic  materials  and  for  such  material  Equations  M2  are  a  valid  ap¬ 
proximation.  Thus  in  all  regions  of  a  mass  filter  that  are  linear,  isotropic,  and 
homogeneous  (i.e.,  within  the  pole  pieces  and  in  the  vacuum  enclosed  by  them)  Maxwell’s 
Equations  M2  are  valid.  In  a  mass  filter,  the  vacuum  conditions  approximate  free  space 
so  that  the  permittivity  e  and  permeability  ft  should  be  well  approximated  by  €(1  and  po. 
Furthermore,  with  a  low  enough  density  of  ions,  p  and  J  will  approach  zero,  and 
therefore  Maxwell's  Equations  M3  are  a  good  approximation  in  the  vacuum  region 
enclosed  by  the  pole  pieces  of  a  mass  filter.  However,  in  a  mass  filter,  the  voltages  on  the 
electrodes  change  at  radiofrequencies  resulting  in  B  and  E  fields  which  will  also  change 
at  the  same  frequencies,  so  it  is  not  obvious  that  Maxwell’s  Equations  M4  are  valid  in  the 
vacuum  region  enclosed  by  the  jkiIc  pieces.  The  use  of  Equations  M4  in  this  region  can  be 
understood  through  the  following  discussion.  Equations  M3  predict”  that  Bu=E(i/c 

'  E.  M.  Pugh  and  F.  W.  Pugh.  Principles  nf  hlectricity  and  Magnetism.  Addison-Wetley  Publishing  Co.,  Reading, 
MA  ( 1960). 
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where  Bn  and  Eo  are  the  maximum  amplitudes  of  the  B  and  E  fields  and  e  is  the  speed  of 
light  where  c  =  (e Jij''*.  Thus  the  B  field  in  the  Lorentz  force  law  (Eq.  2)  exerts  a  negligi¬ 
ble  force  on  the  ion  compared  to  the  E  field  unless  the  ion  is  moving  at  a  speed  com¬ 
parable  to  the  speed  of  light.  The  fastest  ions  in  a  quadrupole  mass  filter  are  moving 
slowly  ( vie  <  10'4)®  compared  to  the  speed  of  light,  even  under  extreme  operating  condi¬ 
tions.  Therefore,  it  is  not  necessary  to  consider  the  B  fields  any  further  for  the  purpose  of 
computing  the  force  on  an  ion.  The  Lorentz  force  law  (Equation  2)  reduces  to 

F  =  eE,  (9) 

while  Maxwell’s  Equations  M4  become: 

V-E  =  0  Vx  E  =  0.  (M5) 

The  second  of  these  Equations  M5  is  valid  only  for  an  alternating  potential  which  has  a 
frequency  v  low  enough  such  that 


X  =  civ  >  >  t  (10) 

That  is,  the  wavelength  X  associated  with  the  electromagnetic  wave  must  be  much  greater 
than  the  length  t  of  the  electrode  structure.  The  necessity  of  Equation  10  can  be 
understood  in  the  following  way.  If  the  wavelength  does  not  satisfy  Equation  10  there  is  a 
possibility  of  standing  waves  (just  as  in  the  case  of  sound  waves  in  an  organ  pipe)  and 
these  are  not  predicted  by  Equations  M5.  Typically,  the  length  of  the  mass  filter  is  less 
than  1  m  and  the  frequencies  used  are  well  below  10  MHz.®  Therefore,  Equations  9  and 
M5  are  valid  for  the  vacuum  region  traversed  by  the  ions  in  the  mass  filter. 


An  E  field  which  satisfies  Maxwell’s  Equations  M5  can  be  found  with  the  aid  of  two 
mathematical  theorems.  Stokes  theorem,  valid  for  any  continuous,  differentiable  vector 
field  A,  asserts  that 


VxA  •  ds 


s 


A  •  dr. 


(ID 


where  dr  is  a  vector  on  the  circumference  of  the  enclosed  area  S  and  whose  direction  is  in 
the  direction  of  the  integration.  The  vector  ds  has  a  magnitude  equal  to  the  area  ds.  and 
its  direction  is  normal  to  ds.  The  ambiguity  in  the  two  directions  normal  to  d>  is  resolved 
by  the  right-hand  rule,  that  is  ds  points  in  the  direction  of  the  thumb  with  the  other 
fingers  of  the  right  hand  pointing  in  the  direction  dr.  In  Equation  1 1.  the  integral  on  the 
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J.  E.  Campana,  Int.  J.  Mass  Spectrom.  Ion  Phys.,  33,  101  (1980)  and  references  cited  within. 


right  is  over  the  entire  eircumferenee  and  the  integral  on  the  left  is  over  the  area  S  en¬ 
closed  by  the  eircumferenee.  Stokes  theorem  enables  a  two-dimensional  integral  over  a 
surface  which  has  the  form  given  on  the  left  side  of  Equation  1 1  to  he  evaluated  by  the 
simpler  one-dimensional  integral  on  the  right-hand  side  of  Equation  11.  Application  of 
Stokes  theorem  to  the  second  of  Equations  M5  implies 

y*E-dr  =  0  (12) 

over  any  closed  path  which  defines  a  surface  over  which  VxE  =  0. 

The  second  mathematical  theorem  asserts  that  anv  vector  field  E.  which  has  the  proper¬ 
ty  that  §E-dr  =  0  around  every  closed  path,  can  he  represented  as  the  divergence  of  a 
scalar  field  <£(x.y.z),  that  is 


E=-V0(x.y.z).  (13) 

The  minus  sign  in  Equation  13  is  not  essential  and  is  included  to  conform  with  conven¬ 
tion.  This  is  analogous  to  the  situation  in  classical  mechanics  where  the  potential  energy 
U(r)  is  defined  so  that  F  =  -VE(r).  The  scalar  field  <£(x.y.z)  is  called  the  potential  func¬ 
tion.  Equation  13  states  that  a  vector  field,  which  normally  needs  three  scalar  fields  to 
represent  it  (one  for  each  component,  can  he  represented  hv  a  single  scalar  field  if  Equa¬ 
tion  12  is  satisfied  for  every  closed  path. 


The  validity  of  Equation  13  in  a  region  of  space  for  which  Equation  12  is  satisfied  for 
every  closed  path  is  seen  by  the  billowing  argument.  Equation  12  applied  to  Figure  1  im¬ 
plies 


y*E-dr  +  yVdr  =  0. 

ABC  CDA 


(14) 


where  the  letters  on  the  bottom  of  the  integral  sy  mbol  indicate  the  path  of  the  integration. 
If  the  direction  of  a  line  integral  is  reversed,  then  the  sign  changes  such  that 


J E-dr  =-/e- 


dr 


CDA  ADC 


and  Equation  14  becomes 


/ 


E  -  dr 


ABC 


ADC 


(15) 
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Figure  1.  Vector  fields  and  gradients. 

If  the  line  integral  of  a  vector  Field  E(x,y,z)  around  any  closed  path  is  zero,  then  there  exists 
a  scalar  quantity  <4(x,y,z)  a  function  of  position  only  with  the  property  E  -  -  V0- 


Equation  15  states  that  if  Equation  12  is  valid  then  the  line  integral  from  A  to  C  is  in¬ 
dependent  of  the  path.  Thus  if  Equation  12  is  valid,  then  the  integral  from  point  A  to 
point  C  can  only  depend  on  the  coordinates  of  points  A  and  C . 


/ 

ABC 


E-dr  =-j«/»(r.)-</»(r4)J. 


(16) 


The  function  <f>  called  the  potential,  corresponds  to  the  voltage  on  the  poles  of  the  elec¬ 
trode  strueture.  Now  suppose  points  A  and  C  are  taken  close  together,  then  Equation  16 
becomes 


E  -  dr  =-d<£  =— V0*< lr,  (17) 

where  the  last  equality  follows  from  the  definition  of  (Eq.  3).  Since  Equation  17  is 
true  for  any  element  of  length  dr,  there  exists  a  function  (/>  which  depends  only  on  the 
postion  such  that  Equation  13  is  true.  This  completes  the  proof  of  Equa¬ 
tion  13. 

Substituting  E  from  Equation  13  into  the  first  equation  of  Equations  M5  gives 

V-V$(x,y.z)  =  0.  (1H) 

The  symbol  V-V  is  called  the  Laplacian  and  is  denoted  by  V2;  consequently.  Equation 
18  becomes 


V2<^(x.y.z)= o.  <i9i 

Equation  19  is  called  Laplace’s  equation  and  in  rectangular  coordinates 

V2  —  V  •  V 

(f ±*i± 

V  dx  9y  dz  /  \  dx  dy  dz/ 

=  il  +  —  *  92 

dx2  dy2  dz2 
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Expressions  for  V2  in  the  cylindrical  and  spherical  coordinate  systems  are  given 
elsewhere.13  14  Since  Equation  19  is  a  consequence  of  Maxwell’s  Equations  M5,  it  is  ap¬ 
parent  that  if  a  function  <f>{\,y,z)  can  be  found  which  satisfies  Equation  19,  then  the  E 
field  defined  by  Equation  13  will  satisfy  Maxwell's  Equations  M5. 

The  prescription  for  finding  the  differential  equation  of  ion  motion  for  any  configura¬ 
tion  of  pole  pieces  can  now  be  given.  Suppose  that  an  ion  is  moving  in  the  field  of  a  set  of 
n  conductors  each  at  a  constant  voltage  (^(i  —  1,2,.. .,n).  l>et  t.  —  f(x.v)  denote  the  equation 
for  the  surface  of  the  ith  conductor.  The  differential  equation  of  motion  can  lie  found  by 
obtaining  a  potential  function  <£(x,y,z)  which  satisfies  Laplace's  equation  (Equation  19) 
and  which  also  satisfies  the  boundary  conditions 

<^(x,y,fi(x,y))=<^i  (i=  1,2... .n).  (20) 

This  latter  requirement  says  that  the  potential  function  must  be  equal  to  the  voltage  ap¬ 
plied  to  any  of  the  electrodes  at  the  electrode  surfaces. 

Generally,  it  is  difficult  to  find  such  a  function  which  satisfies  Equations  19  and  20. 
The  function  0(x,y,z)  can  sometimes  be  found  analytically  using  the  method  of 
images,13  14  17  complex  variable  theory.13  14  17  13  infinite  series  expansions13  14  17  14  or  it 
can  be  evaluated  numerically  using  digital  computers.20  21 


E.  M.  Pugh  and  E.  W.  Pugh,  Principles  of  Electricity  and  Magnetism,  Addison-Wcsley  Publishing  Co.,  Reading. 
MA(1 960). 
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S.  Ramo,  J.  R  Whinncry,  and  T.  Van  Duzer,  Fields  and  Waves  in  Communication  Electronics.  John  Wiley  and 
Sons,  Inc.,  New  York  (1965). 

*  J.  D.  lackson.  Classical  Electrodynamics,  John  Wiley  and  Sons,  Inc.,  New  York  (1962). 

18  R.  V.  Churchill,  Complex  Variables  and  Applications,  McGraw-HiJJ  Book  Co.,  Inc.,  New  York  <  I960). 

R.  V.  Churchill,  Fourier  Series  and  Boundary  Value  Problems.  McGraw-Hill  Book  Co.,  Inc.,  New  York  (1941). 

^  K.  J  Binns  and  P.  J.  Lawrenson,  Analysis  and  Computation  of  Electric  and  Magnetic  Field  Problems.  The 
Macmillan  Co.,  New  York  (1963). 

M.  L.  Bar on  and  M.  G.  Salvador),  Numerical  Methods  in  Engineering,  Prentice-Hall,  Inc.,  Englewood  Cliffs, 
NJ  (1964). 
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Once  the  potential  dKx.v.z)  is  known,  the  differential  equation  of  motion  can  lie  found 
using  Newton's  law  of  motion 


m  =  -  c  V#x,  y.  z). 

dtJ 


(21) 


As  discussed  earlier,  this  method  is  approximately  valid  for  time-varying  potentials** 


=  '‘os  W,‘ 


(22) 


providing  the  potential  does  not  change  too  rapidly  with  time  (i.e.,  subject  to  the  limita¬ 
tions  of  Equations  M3  and  10).  Physically.  Equation  22  corresponds  to  subjecting  the 
electrodes  to  a  voltage  which  has  both  d.c.  and  a.c.  components. 

The  remainder  of  this  report  will  consider  only  infinitely  long  conductors  whose  shape 
d<tes  not  vary  with  z.  Although  real  electrode  structures  must  necessarily  be  of  finite 
length,  this  approximation  is  made  because  it  captures  the  essential  features  for  the 
o|ieration  of  electric  RF  devices  while  avoiding  some  of  the  mathematical  complexities. 
With  this  approximation,  the  potential  for  such  a  conductor  configuration  does  not  de¬ 
pend  on  z;  i.e..  (tftx.v.z)  =  0<x.v)  and  hence  the  force  in  the  z-direction  is  zero;  i.e.. 


dJz 


dtJ 


=  0.  The  equation  of  motion  in  the  z-direction  can  be  integrated  directly 


z  =  v  t  +  z 

t  <• 


(23) 


Here  v  ami  zm  are  the  z-component  of  velocity  and  the  z-coordinate  respectively  at  t  =  (). 

It 

III.  METHOD  FOR  OBTAINING  POTENTIAL  FUNCTIONS  AND 
EQUATIONS  OF  ION  MOTION  IN  MULTIPOLE  FIELDS 

The  approach  of  the  previous  section  was  to  find  the  potential  function  dHx.v.z)  for  a 
sfiecified  set  of  conductors  at  potentials  <t>r  An  alternative  approach  to  the  problem  is  to 
find  a  function  IJx.y)  which  satisfies  Laplace's  equation  (Equation  19)  and  from  this 
function  determine  the  equipotential  surfaces  of  the  pole  pieces.  If  conductors  with  given 
applied  jiotentials  are  fashioned  which  coincide  with  the  equipotential  surfaces,  a  poten¬ 
tial  function 

^  The  minus  sign  in  hquation  22  is  not  essential  and  is  chosen  by  convention. 


^(x.y|=KIiIJn(x,y)  (24) 

<*an  lie  found  which  satisfies  Laplace's  equation  and  the  boundary  conditions  (Kquation 
20)  hy  simply  choosing  the  constant  Kn  properly.  Physically.  Kn  is  related  to  the  dimen¬ 
sion  of  the  electrode  structure  and  the  applied  electrode  potentials;  it  is  chosen  such  that 
$(x,y)  matches  the  known  potentials  at  the  equipotential  surfaces. 

Functions  which  satisfy  Laplace’s  equation  may  be  found  by  application  of  the  theory 
of  complex  variables.  If  the  complex  variable  z  =  x  +  iy  is  raised  to  an  integral  power  n. 
then  the  result  can  lie  expressed  as  the  sum  of  two  functions,  a  real  Un(x.v)  and  an 
imaginary  Vn(x,y)  part 


(x  +  iy)"  =  Un(x,v>  -I-  iVn(x,y).  (25) 

The  integer  n  defines  the  order  of  the  multipole  field.  From  complex  variable  theory  z"  is 
known  to  lie  analytic  and  so  the  (lauchy-Riemann  equations  are  applicable:18 


dUn(x,y)  _  aVn(x,y) 
dx  dy 


3Un(x,y)  dVn(x,y) 

dy  dx 


(26a) 


(26l>) 


The  sum  of  the  partial  derivative  of  the  first  of  these  equations  with  respect  to  x  and  the 
partial  derivative  of  the  second  of  these  equations  with  resjiect  to  y  gives 


d1  Un(x.y) 
dx2 


3J  Un(x,y) 

d? 


=  0. 


(27) 


Thus  the  function  lln(x,y|  satisfies  Laplace's  equation. 

The  method  descrilied  here  for  obtaining  the  potential  function  and  the  differential 
equation  of  motion  can  Ite  summarized  as  follows. 

R.  V  Churchill.  Complex  Variables  and  Applications,  McCuaw-Hill  Book  Co  ,  Inc..  New  York  ( I960). 
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A.  Evaluate  the  funetion  Un(x,y)  for  the  variable  z“  where  n  is  any  positive  non-zero  in¬ 
teger.  The  funetion  Un(x,y)  will  satisfy  Laplace’s  equation  (Equation  27). 

B.  The  equipotential  conductor  surfaces  are  constructed  so  that  they  fall  on  the  locus 
of  points  in  the  (x,y)  plane  defined  by  Un(x,y)=  ±  constant.  The  separation  between  op¬ 
posite  electroties  is  chosen  to  be  2r  ( to  insure  that,  for  the  case  of  the  quadrupole,  the  stan¬ 
dard  quadnipole  geometry2  is  reproduced.  Geometrically  adjacent  conductors  are  ar¬ 
bitrarily  chosen  to  have  applied  potentials  of  0.  =  <f>J 2  and  4>i  =  -0<>/2,  where  0o  is 
allowed  to  vary  with  time  (Equation  22),  subject  to  those  constraints  put  forth  previously 
(the  constraints  of  Equation  10  and  the  constraints  associated  with  Equations  M5).  With 
the  0.  so  defined,  0  is  simply  the  voltage  between  any  two  adjacent  electrodes. 

(1.  An  appropriate  constant  K  is  determined  such  that  KnUn(x,y)  =  0(x,y)  satisfies 
the  boundary  conditions;  i.e.,  Kn  is  chosen  so  that  0(x,y)  matches  the  known  potentials  at 
the  equipotential  surfaces. 

I).  The  differential  equation  of  an  ion  is  determined  from  the  potential  found 
in  (]  using  Newton’s  law  of  motion  (Equation  21). 

In  brief,  find  a  potential  function  0(x, y)  which  satisfies  Laplace’s  equation  (Equation  19) 
and  the  boundary  conditions  (Equation  20);  then  ion  motion  is  found  from  Newton's  law 
of  motion  (Equation  21). 

This  method  will  be  applied  in  Section  V  to  obtain  potential  functions  and  differential 
equations  of  motion  for  various  muliipole  electrode  geometries. 

IV.  APPROXIMATIONS  IN  THE  EQUATIONS  FOR  ION  MOTION 

Several  approximations  implicit  in  the  prescription  for  finding  the  differential  equa¬ 
tions  of  ion  motion  given  in  the  preceding  section  are  a  result  of  neglecting  the  various 
forces  that  ions  may  experience.  These  “idealized”  equations  of  motion  are  normailv  used 
in  describing  ion  motion.2  7'10  These  approximations  will  now  be  explicitly  discussed. 

^  W.  Paul  and  H.  Steinwcdcl,  Naturforsch.  Teil  A.,  8.  448  (1953);  W.  Paul,  H.  P.  Reinhard  and  von  Zalin. 

Z.  Phys.  152,  43(1958). 

7  P.  11.  Dawson  (Editor),  Quadrupole  Mass  Spectrometry  and  Its  Applications,  Elsevier.  Amsterdam  (1976). 

fi 

J,  I  J.  Todd  and  G.  Lawson.  MTP  International  Review  of  Science,  Physical  Chemistry,  Mass  Spectrometry. 

Series  Two,  Vol.  5,  edited  by  A.  Maccoll,  Buttcrworths,  London  (1975). 

o 

J.  F..  C’ampana,  Int  J.  Mass  Spectrom.  Ion  Phvs.,  dd,  101  (1980)  and  references  cited  within. 

' I).  R  Denison,/  Vac.  Sci.  Technol..  8,  266  (1971). 
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A.  Finite  Length  of  the  Electrode  Structure.  The  differential  equations  were 
derived  under  the  assumption  that  the  electrode  structure  is  infinitely  long.  Real  electrode 
structures  are  not  infinitely  long  and  so  the  electromagnetic  fields  are  not  functions  of  x 
and  y  alone  hut  also  depend  on  z;  i.e.,  they  do  not  begin  and  end  abruptly  at  the  ex¬ 
tremities  of  the  electrode  structure  but  extend  at  either  end  of  the  electrode  structure  giv¬ 
ing  rise  to  what  is  termed  fringing  fields.  Fringing  fields  have  been  neglected  in  these 
derivations  but  for  an  actual  mass  filter  these  fields  may  be  important.7  The  assumption 
of  infinite  length  was  implicitly  made  when  the  potential  was  assumed  to  be  independent 
of  the  z-coordinate. 

B.  Image  Force  on  ions.  If  an  ion  with  a  charge  e  is  placed  near  an  infinite 
conducting  she'-t  U  ground  potential,  the  field  of  the  ion  will  cause  the  free  charges 
within  the  conoV'tor  to  move  in  such  a  way  that  the  field  within  the  conductor  is  zero. 
These  induced  n >  on  the  conductor  will  attract  the  ion  to  the  conductor  with  a 
force,13 


where  d  is  the  distance  from  the  charge  to  the  conducting  sheet.  It  is  as  if  there  were  an 
image  charge  of  opposite  polarity  on  the  other  side  of  the  conductor,  thus  the  name 
“image  force”  is  used.  Similarly,  a  charge  in  a  conducting  hollow  sphere  (not  at  the 
center)  will  he  attracted  to  the  nearest  point  of  the  conductor.13  For  the  geometry  of  the  in¬ 
finite  sheet  or  sphere,  the  force  on  the  ion  is  non-zero  and  this  suggests  that  the  force  for 
an  ion  moving  in  the  multipole  geometries  considered  here  is  non-zero.  The  forces  due  to 
the  induced  charges  on  the  conductor  have  been  neglected  in  going  from  Equations  M2  to 
M3:  the  radiofrequency  source  allows  charge  to  flow  to  and  from  the  electrode  structure 
and  so  Equation  M3  is  not  valid  in  the  region  of  space  which  includes  the  electrode  struc¬ 
ture.  In  typical  computer  trajectory  calculations,  the  ion  is  assumed  to  be  lost  once  it  hits 
the  electrode.23  Because  of  this  so  called  “image  force,”  the  ions  are  more  properly  lost 
when  they  get  to  within  a  certain  distance  of  the  electrode  structure. 

E  Fields  Induced  by  Ions.  Electric  forces  between  ions  may  occur  within 
the  vacuum  region  enclosed  by  the  pole  pieces  when  an  ion  creates  an  electric  field  E  at 
the  site  of  a  second  ion.  The  second  ion  is  thus  acted  upon  by  an  electric  force  according 
to  Coulomb's  law. 

7  P.  H.  Dawson  (Editor),  Quadrupolc  Mass  Spectrometry  and  Its  Applications,  Elsevier,  Amsterdam  (1976). 

E.  M,  Pugh  and  F..  W.  Pugh,  Principles  of  Electricity  and  Magnetism,  Addison-Wesley  Publishing  Co.,  Reading. 

MA  (I960). 

^  J.  E.  Campana  and  P.  C.  Juts,  Int.  J.  Mass  Spcctrom.  Ion  Phys.,  .13,  1 19  (1980). 


D.  B  Fields  Induced  by  Ions.  The  movement  of  one  ion  generates  a  B  field 
associated  with  the  movement  of  that  ion  and  if  a  second  ion  moves  in  the  B  field 
generated  by  the  first  ion,  it  will  feel  a  force  due  to  the  B  field  of  the  first  ion  providing  v 
is  not  parallel  to  B.  This  force  is  neglected  in  the  derivation  given  here.  The  forces  in  B, 
C,  and  D  were  neglected  when  J  and  P  were  set  equal  to  zero  in  Section  I  in  going  from 
Equation  M2  to  M3.  These  latter  two  forces  discussed  in  C  and  D  can  be  collectively 
termed  space  charge  effects. 


E.  Presence  of  B  Fields  on  Computing  the  E  Field.  The  E  field  computed  from 


the  potential  is  only  approximately  correct  because  the  term  with _ has  been  negleeted 

dt 


in  the  second  equation  of  Equation  M5.  At  high  enough  frequencies,  or  with  a  mass  filter 
with  very  large  dimensions,  the  E  field  computed  from  the  potential  will  be  significantly 
in  error  and  this  would  result  in  significant  error  in  the  differential  equation  describing 
the  ion  motion. 


F.  Presence  of  B  Field  on  Computing  the  Force  of  the  Ion.  The  fourth  equation 
of  Equations  M3  shows  that  a  changing  E  field  generates  a  B  field  and  from  the  Lorentz 
force  law  this  B  field  exerts  a  force  on  the  ion.  Unless  the  ion  is  moving  close  to  the  speed 
of  light,  this  force  is  usually  negligible  compared  to  the  E  field.  It  is  reasonable  to  neglect 
this  term  for  practical  applications. 

C.  Radiation  by  the  Accelerated  Ion.  A  consequence  of  Maxwell's  equations  is 
that  an  accelerating  charge  gives  off  electromagnetic  radiation.17  An  accelerating  charge 
requires  an  effective  force  to  help  balance  the  energy  radiated.  This  force,  which  is  pro¬ 
portional  to  the  acceleration  of  the  ion,17  has  been  neglected  in  these  derivations. 

The  approximations  discussed  in  E,  F,  and  G  were  introduced  as  a  consequence 
of  dropping  the  time  dependence  of  the  E  and  B  fields  in  going  from  Equations  M3  to 
M4. 


H.  Ion  Neutral  Interaction.  In  the  normal  operation  of  a  mass  filter,  ions  may  in¬ 
teract  with  residual  gas  in  the  chamber.  It  has  been  suggested24  that  this  force  can  be 
represented  as  a  viscous  drag  on  the  motion  of  the  ion.  The  viscous  force,  proportional  to 
the  speed  of  the  ion  has  been  neglected  in  the  derivation  of  the  equations  given  here.  At 
the  normal  operating  pressures  of  mass  filters  (10's  mm  Hg).  the  mean  free  path  of  a  par- 
tide  is  greater  than  4.0  x  10'3  m. 

17  J.  D.  Jackson,  Classical  Electrodynamics,  John  WiJey  and  Sons,  Inc.,  New  York  (1962). 

N.  R.  Whetten,  J.  Vac.  Sci.  Technol.,  II,  515  (1974). 


16 


I.  Gravitational  Interaction.  The  gravitational  interaction  has  been  neglected. 
Because  of  the  small  amount  of  time  an  ion  spends  traversing  a  quadrupole  mass  filter, 
this  effect  is  small  for  such  a  device  and  is  expected  to  he  more  important  in  ion  storage 
devices  (three-dimensional  mass  filter).7 

J.  Relativistic  Effects.  Non-relativistic  laws  of  mechanics  are  used  in  the  method 
described  in  Section  II.  Since  the  fastest  ions  in  a  quadrupole  mass  filter  move  slowly9 
(v/c  <  10"4)  compared  to  the  speed  of  light,  the  errors  introduced  by  using  non-relativistic 
mechanics  are  expected  to  be  small. 

V.  APPLICATION  TO  SOME  MULTIPOLE  GEOMETRIES 

The  method  described  in  Section  III  will  be  demonstrated  here  to  obtain  the  potential 
functions  and  the  differential  equations  of  ion  motion  for  the  dipole,  quadrupole,  hex- 
apole,  octapole,  and  decapole  electrode  geometries. 

A.  The  Dipole  Field.  The  variable  z  is  raised  to  the  first  power  to  obtain: 

U,(x,y)  =  x 
vi(x.y)  =  y. 

SinceV2lJ1(x,y)  =  0.  the  function  U,(x.y)  forms  the  basis  for  a  potential  function  in  charge- 
free  space.  Its  equipotential  surfaces  are  illustrated  in  f  igure  2.  This  leads  to  a  model  in 
which  there  are  two  plane  sheets:  one  at  x=r  and  one  at  x  =  -r  with  applied  potentials  of 
4>J2  and  -<jW2.  respectively.  The  potential  function  0(x.v )  =  K,L',(x.y)  can  be  deduced  to 
be 

X 

0(x.y)  = - :  -ro  <x  <ro 

2  r 


if  it  is  to  satistv  both  Laplace's  equation  and  the  boundary  conditions 
|</>(  ±  r (  .y)=  db  <f>  /2|  in  Equation  20.  With  given  by  Eq.  22.  Newton's  law  Eq.  21 
becomes: 
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Figure  2.  Dipole  geometry. 

The  positive  and  negative  equipotential  surfaces  consists  of  two  plane  sheets  separated 
by  a  distance  2ro.  The  surfaces  at  x»r0  and  x-  -  rD  are  at  potentials  <p0/ 2  and  -  <t>a/2. 


i. 

s 


18 


|  nfi  WMi 


B.  The  Quadrupole  Field.  The  function  z"  for  n  =  2  is  evaluated  to  obtain: 

U2(x,v)  =  x2  —  y2 
V2(x.y)  =  2xy. 

SinceV2U2(x,y)  =  0,  the  function  U2(x,y)  is  the  basis  for  a  possible  potential  function  in 
charge-free  space.  Its  equipotential  surfaces  are  illustrated  in  Figure  3.  Each  of  the  four 
hyperbolic  pole  pieces  have  the  same  geometry  and  hence  the  structure  is  unchanged 
under  90°  rotation.  The  pole  pieces  lying  on  the  x-axis  are  arbitrarily  chosen  to  have 
positive  applied  potentials  +<j>J2  and  the  pole  pieces  lying  on  the  v-axis  are  arbitrarily 
chosen  to  have  negative  applied  potentials -0  /2.  The  potential  function  0(x.v)  can  he 
deduced  to  be 


0(x,y)  =  _ 1_  |x2  —  v2). 

Jr  2 


It  is  readily  verified  that  this  potential  function  satisfies  both  Laplace's  equation  and  the 
boundary  conditions.  Equation  20  (eg.  0(r  .0)  =  0  /2  and  0(O.r  =  —  0  /2].  With  0 
given  by  Equation  22,  Newton’s  law  becomes: 


d2  x 

m _  + 

dt2 


(0[)C  -  0AC  cos  cot)  x  =  0 


m -  -  -  (0^  -  0A(,  cos  cot)  y  =  0 


The  electric  potential  and  the  differential  equation  obtained  here  agree  exactly  with  those 
derived  elsewhere  by  various  other  methods.*  Equations  29  are  of  the  well  studied 
Mathieu-tvpe  equation.4"6 

4 

f'.  M.  Arscott,  Periodic  Differential  Equations.  The  Macmillan  C  o..  New  York  ( l%4). 

^  N.  W.  McLachlan,  Theory  and  Application  of  Mathieu  Functions.  Oxford  University  Press.  New  York  11947). 

^  Computation  Laboratory.  Tables  Relating  to  Mathieu  Functions.  C.S.  Bureau  of  Standards.  Columbia  Univer¬ 
sity  Press,  New  Y  ork  (1951). 
q 

J.  I*.  Campana,  Jnt.  J.  Mass  Spcctrom.  Ion  Phvs.,  33.  101  ( 1980)  and  references  cited  within 
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Figure  3.  Quadrupole  geometry. 

The  positive  and  negative  equipotentials  consist  of  four  surfaces  with  an  inscribed 
radius  of  rQ.  The  equations  of  the  equipotential  surfaces  are  0(x,y)=(0o/2ro:  )(x‘-y: ). 
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('.  The  Hexapole  Field.  The  fund  ion  z"  for  n  =  3  is  evaluated  to  obtain: 


L:.,(x.y)  =  xs  —  3xy2 
Vt(x.v)=  3x*v  —  y\ 


SineeV?U.t(x.y>=  0,  the  funetion  U.,( x.y )  is  the  basis  for  a  possible  potential  function  in 
charge-free  space.  Its  equipotential  surfaces  are  illustrated  in  Figure  4.  Each  of  the  six 
pole  pieces  have  the  same  shape  and  so  the  structure  is  unchanged  under  a  60°  rotation. 
Similarly,  as  in  the  previous  examples,  the  potential  function  y)  can  be  deduced  to  be 


<j> 

4>(x,y)  =  — —  (x3  -  3xy2 ). 
2r  3 


It  is  readilv  verified  that  this  function  satisfies  both  Laplace's  equation  and  the  boundary 
conditions.  With  </>u  given  by  Equation  22.  Newton's  law  becomes: 


m . 


d2x  +  3  e 
dt2  2  r  - 


(<t>0 c  -  0Ac  cos  (x2  -  y2 )  =  o 


d2y  3e  ,,  _ 

m _ i  -  _  ((Pqc  -  <pAC  cos  tot)  xy  =  0. 

dt2  r  3 


(30) 


D.  The  Octapole  Field.  The  function  z"  for  n  =  4  is  evaluated  to  obtain: 

Ut(x.y)=  x4  —  6x2v2  +  v* 

Vj(x.v)=  4xy(x2  —  y2). 

Since V2l!f(x,y)  =  0,  the  function  U((x.y)  is  the  basis  for  a  possible  potential  function  in 
charge-free  space.  Its  equipotential  surfaces  are  illustrated  in  Figure  5.  Each  of  the  eight 
pole  pieces  have  the  same  shape  and  so  the  structure  is  unchanged  under  a  45°  rotation. 
Similarly,  as  in  the  previous  examples,  the  potential  d>(x,y)  can  be  deduced  to  be 


0, 


0(x,y)  = _  (x4  -  6x2y2  +  y4  ). 

2r  4 
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Figure  5.  Octapole  geometry. 

The  positive  and  negative  equipotentials  consists  of  eight  surfaces 
with  an  inscribed  radius  of  rQ.  The  equation  of  the  equipotential 
surfaces  are  tf>ix,y)  -  (0o/2ro4)  (x4  -  6x3y2  +  y4). 
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It  i.s  readily  verified  that  this  function  satisfies  Laplace's  equation  and  the  boundary  con¬ 
ditions.  With  <f>t)  by  Equation  22,  Newton’s  law  becomes: 

m  +  ~e.  (0^  -  0AC  cos  cot)  (x3  -  3xy2 )  =  0 
dt2  r  4 

O 

.2  -y  (31) 

m  — X  -  __L  (0^  -  0AC  cos  cot)  (3x2  y  -  y3 )  =  0. 
dt2  ro4 

E.  The  Decapole  Field.  The  function  z"  for  n  =  5  is  evaluated  to  obtain: 

U_(x,y)=xs  —  10x3ys  +  5xv4 
V5(x,y)=ys  —  10x2y3  4-  5x4y. 

Since  V2U.(x,y)  =  0,  the  function  U5(x,y)  is  the  basis  for  a  possible  potential  function  in 
charge-free  space.  Its  equipotential  surfaces  are  illustrated  in  Figure  6.  Each  of  the  10 
pole  pieces  have  the  same  shape,  and  so  the  structure  is  unchanged  under  a  36°  rotation. 
Similarly,  as  in  the  previous  examples,  the  potential  function  <0(x,y)  can  be  deduced  to  be 

0 

0(x,y)=  — (xs  -  10x3y2  +  5xy4). 

V 

It  is  readily  verified  that  this  function  satisfies  Laplace’s  equation  and  the  boundary  con¬ 
ditions.  With  <0  given  by  Equation  22,  Newton’s  law  becomes: 

m  +  — — -  (0dc  -  0Ac  cos  wt)  (x4  -  6x2y2  +  y4 )  =  0 
dt2  2ro  5 

(32) 

in  (0™-  -  0. „  cos  coy)  (x3y  -  xy3 )  =  0. 

dt2  r  5 

O 

F.  Fields  of  Geometry  Higher  Than  the  Decapole.  The  method  illustrated  in  these 
examples  can  be  used  without  difficulty  for  finding  the  differential  equation  of  an  ion 
for  n  >  3.  where  n  is  an  integer. 

G.  Generalization  of  Differential  Equations  to  an  Arbitrary  Phase.  The  alter¬ 
nating  component  at  time  t  =  0  is  0A(.  in  the  differential  Equations  28-32.  Ions  may  enter 
the  electric  field  at  any  time  where  the  alternating  component  can  have  any  value 
between  +  <0A(  and-<0A) .  Thus  Equations  28-32  may  be  written  more  generally  by  replac¬ 
ing  cat  with  cat  +  6  where  ()<5< 2ir. 


X 

Figure  6.  Decapole  geometry. 

The  positive  and  negative  equipotentials  consist  of  10  surfaces  with  an  inscribed  radius 
of  rn.  The  equations  of  the  equipotential  surfaces  are  <Mx,y)=(0o/2ro5  )(xs-10x3yJ+5xy4 ). 


25 


H.  Redundancy  of  the  V  Function.  The  function  V  (X.v)  defined  in  Kquation  25 
also  satisfies  Laplace's  equation:  i.e..V2Vn(x.v  )  =  0.  This  can  he  seen  readily  h>  taking  the 
sum  of  the  partial  derivative  of  Kquation  26a  with  respect  to  \  and  the  partial  derivative 
of  Kquation  26h  with  respect  to  x.  Since  \m(x.v>  satisfies  Laplace's  equation,  it  might 
seem  that  the  equipotentials  Vjx.y)  =  <t>J2  and  YJx.v)  =-0/2  would  generate  different 
surfaces  than  the  corresponding  equations  for  the  function  l  n(x.y). 

The  surfaces  defined  by  Vjx.v)  =  ±  0„/2  are  the  same  as  those  defined  by 
l  Jx.v)  =  ±  0/2  respectively  except  rotated  by  an  angle  a  .  This  is  seen  bv  expressing 
(x  +  iy|"  in  polar  coordinates 

(x  +  iy>"  =  (re'V  =  r"(cos  nd  +  i  sin  nd). 

( .omparison  of  this  last  equation  with  Kquaiton  25  shows 

lJ„(x.v)  =  rn  cos  n0  (33a) 

V„(x,y)=rn  sin  110.  (33l») 

Define  a  primed  coordinate  system  rotated  by  an  amount  an=ir/2n  relative  to  the  un- 

primed  O'  =  0  -  — 

In 


In  the  primed  coordinate  system 

Vn  =  rn  sin  n (O'  +  JL 

\ 

Since  the  form  of  Kquation  54  is  the  same  as  Kipiation  55a.  the  use  of  the  function 
Vn(x.y>  generates  the  sane-  apparatus  as  the  function  l  nlx.v|  rotated  bv  an  angle  at  =  n 

2n 

VI.  Ml  LTIPOLK  <;K()MKTRIKS  AS  MASS  HI.TKRS 

The  methods  presented  in  the  preceding  sections  yield  the  differential  equations  of  mo¬ 
tion  for  ions  in  time-vary  ing  electric  fields.  A  system  functions  as  a  mass  filler  if  it  passes 
ions  in  a  certain  mass  interval,  rejects  all  other  ions,  and  the  mass  interval  is  under  the 
control  of  the  experimenter. 


r  cos  n O' 


1 54 1 
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A.  The  Dipole  Field.  The  equation  of  motion  for  the  x  and  y  directions  can  be  in* 
Unrated  directly: 


tip rv,  ei.r  COS  U)t 

X  =  A+Bt-  tJ  -  AC 

4mr  2mr  us 

O  O 

(35) 

y  =  vy  t  +  y0- 

o 

Here  A  and  B  are  constants  which  depend  on  the  initial  x-coordinate  and  x-velocitv. 
respectively,  and  yo  and  are  the  initial  v-coordinate  and  y-velocity.  It  is  apparent  from 

Equations  35  that  the  motion  is  unstable  in  both  the  x-  and  y-coordinate  directions  for  all 
initial  conditions.  Thus,  the  dipole  field  with  sinusoidal  voltage  is  unsuitable  as  a  mass 
filter. 


B.  The  Quadrupole  Field.  It  is  well  known  that  once  can  choose  <p[M:  and  <PM  such 
that  there  is  stability  for  both  the  x-  and  ydirections  for  a  narrow  range  of  mass  values. 
The  quadrupole  geometry  is  suitable  as  a  mass  filter.7'9 

C.  The  Hexapole,  Octapole.  and  Decapole  Fields.  The  differential  equations  for 
these  three  geometries  are  not  readily  integrated.  It  has  been  shown  that  numerical 
methods  can  be  used  successfully  to  integrate  equations  of  motion,  when  motion  in  one 
coordinate  direction  is  not  independent  of  the  motion  in  the  other  direction. 10  The  utility 
of  such  systems  as  mass  filters  awaits  the  numerical  integration  of  the  equation  of  motion 
and  is  the  subject  of  future  research. 

VII.  SUMMARY 

A  general  method  for  finding  the  differential  equations  of  ion  motion  in  elec¬ 
tromagnetic  fields  has  been  presented  and  applied  to  five  different  multipole  geometries. 
The  fundamental  physics  and  the  approximations  used  in  obtaining  the  differential 
equations  have  been  reviewed  so  that  limitations  in  the  analytical  theory  can  be  ap¬ 
preciated  better. 

7  P.  H.  Dawson  (Editor),  Quadrupole  Mass  Spectrometry  and  Its  Applications.  Elsevier.  Amsterdam  ( 1 976). 

^  J.  E.  J.  Todd  and  G.  Lawson,  MTP  International  Review  of  Science.  Physical  Chemistry.  Mass  Spectrometry, 
Series  Two,  Vol.  5,  edited  by  A.  Maccoll,  Butterworths,  London  (1975). 

^  J.  E.  Campana,  Ini.  J.  Mass  Spectrom.  Ion  Phys..  33.  101  ( 1980)  and  references  cited  within. 

I).  R.  Denison,/  Vac.  Sci.  Technol.  8.  266  (1971). 
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Today  the  quadrupole  mass  analyzer  is  the  most  widely  used  mass  analyzer  for  low 
resolution  applications.  Typically,  mass  filters  are  used  in  laboratory  environments  where 
simplicity,  compactness,  economy,  absence  of  magnetic  fields,  and/or  the  capability  of 
fast  scan  rates  (especially  for  chromatographic  combinations)  are  important.  Thus,  an  im¬ 
provement  in  resolution/transmission  qualities  of  electric  RF  mass  filters  would  enable 
this  widely  used  mass  spectrometer  to  take  a  larger  role  in  analytical  laboratories  and  in 
basic  research. 

Other  design  considerations  besides  resolution/transmission  characteristics  are  impor¬ 
tant.  If  mass  filters  could  be  made  which  are  smaller,  more  rugged,  and  less  susceptible  to 
vibration  with  reduced  power  requirements  while  being  maintained  more  easily,  thev 
might  find  application  outside  the  laboratory.  Possible  applications  include  surveying  for 
natural  resources,  pin-pointing  sources  of  industrial  and  natural  pollutants  and  the  detec¬ 
tion  of  explosive  materials. 
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ABBREVIATIONS  AND  SYMBOLS 


A 

B 

B 

V 

D 

<1 

E 

E, 

e 

F 

H 

i 

J 

A 

J 

K 

t 

( 

M 

m 


Any  vector  field 
Magnetic  field  vector 
Maximum  amplitude  of  the  B  field 
Speed  of  light 

Electric  displacement  vector 

Distance  from  a  point  charge  to  a  conducting  sheet 

Electric  field  vector 

Maximum  amplitude  of  E  field 

Electron  charge 

Force  vector 

Magnetic  field  intensity  or  magnetic  field  strength  vector 
l  nit  vector  along  x-axis 
Current  density  vector 
l  nit  vector  along  y-axis 

A  constant,  related  to  the  dimensions  of  the  electrode  structure  and  Un¬ 
applied  electrode  voltage 

l  nit  vector  along  z-axis 

length  of  the  electrode  structure 

Magnetization  vector 

Magnetic  polarization  vector 
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t 


ni 


Mass 


Polarization  vector 


Radiofrequency 


Radius  vector 


Inscribed  radius  of  the  electrode  structure 


A  real  function  satisfying  Laplace's  equation 
A  potential  energy  function 

An  imaginary  part  of  a  function  that  satisfies  Laplace's  equation 
Velocity  vector 


Radian  measure 


Radian  measure 


The  permittivity 

The  permittivity  of  free  space 


Radian  measure 


Wavelength 

Permeability 

Permeability  of  free  space 
Frequency 


4Hx.v.z) 


Charge  density 
A  potential  function 


4>  A  scalar  function  of  position  or  applied  voltage 

X,  Polarizability  or  the  electric  susceptibility 

Xn  Magnetic  susceptibility 

U)  Angular  frequency 

V  “del  operator” 

V2  The  Laplacian 

d  Backcurling  delta  (partial  derivative  sign) 

$  Contour  integral 
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